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Abstract 

We study the non-linear Schrodinger equation in (1 + 1) dimensions in which the nonlinear term is 
taken in the form of a nonlocal interaction of the Coulomb or Yukawa-type. We solve the equation 
numerically and find that, for all values of the nonlocal coupling constant, and in all cases, the 
equation possesses solitonic solutions. We show that our results, for the dependence of the height 
of the soliton on the coupling constant, are in good agreement with the predictions based on an 
analytic treatment in which the soliton is approximated by a gaussian. 
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I. INTRODUCTION 



The Nonlinear Schrodinger equation (NLS) in (1 + 1) dimensions arises in many descrip- 
tions of physical phenomena which range from the water waves and wave propagation in 
fibre optics to biophysics Moreover, being perhaps the simplest nonlinear PDE, it has 
been studied a lot by both mathematicians and physicists. In particular the equation was 
shown to be integrable and that it has soliton-like solutions jH . Most of this work involved 
the simplest version of the equation - in which the nonlinearity was in the form of a simple 
power of the field itself. 

The equation has also been extended in many ways; people looked at the effects of the 
additions of higher powers of the fields and of their derivatives. Some of these generalisations 
involved also nonlocal terms. In fact, some of them, in the optics literature, go under the 
name of the Gross-Pitaevskii equation 0|. Due to its general nature the equation has been 
used in many applications. Recently these generalisations involving a nonlocal nonlinear 
term have been used in the description of the Bose-Einstein condensates below the transition 
temperature In this case the external potential models the trap for the particles, while 
the nonlocal interaction describes the interaction between the bosons in the condensate. The 
nonlocal equation has also recently been used to describe wave dynamics in media with a 
large response length, e.g. in nonlinear optics (see -5] and references therein). The nonlocal 
equation was also studied some time ago by Ginibre and Velo who, in their mathematical 
paper, proved that some nonlocal Schrodinger equations possess well defined solutions. This 
proof did not give any details of their nature, i. e. it did not answer the question whether 
the solutions represent solitons and what their properties are. 

More recently Abe and Ogawa looked at the nonlocal version of the nonlinear 
Schrodinger equation (1 + 1) dimensions for a class of nonlinear terms involving a power-law 
interaction. They demonstrated the existence of solitonic solutions and then, using a gaus- 
sian approximation, studied the dependence of the width of the solitons on the exponent 
of the power-law interaction. However, the study did not include the cases of a Coulomb 
or Yukawa interactions. In [?| the study was extended to more than one spatial dimension 
using an exponential-type nonlocal interaction. 

Given that in most applications the nonlocal term comes from the interaction of some 
(spread-out) physical charges we have decided to have another look at this equation for the 
cases of Coulomb - and Yukawa-like interactions. 

In this paper, we present the results of our studies of a system in one spatial dimension 
with nonlocal interaction, i.e. when the kernel of the nonlocal term is given by — k\x — y\ or 
^, exp(-7|x-y|) ^ p a p er j g organised as follows: In Section II, we introduce all the notations 
and present the model. In Sections III and IV, respectively, we discuss the cases of the 
Coulomb - and Yukawa interactions. In each Section we present the numerically obtained 
results for the existence of solitons and their properties and then present some analytic 
arguments based a gaussian approximation which allow us to explain these results. We 
finish with some conclusions in Section V. 
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II. THE MODEL 



The action of the model in (1 + 1) dimensions, which we have studied, is given by: 
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and the corresponding equation reads: 
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Here K(x, x') is the kernel of the non-local interaction of the system. For K(x, x') = 5(x—x'), 
i.e. in the "local" limit, (j2J) reduces to the usual non-linear Schrodinger equation, which 
is known to possess soliton solutions, while the case K(x,x') — \x — x'\~ k , < k < 1 
was studied in p. Here, we are interested in the physically more relevant Coulomb and 
Yukawa- type interactions which lie outside the range of potentials studied in 

Note that the system possesses some conserved quantities, like the energy which is given 

by: 
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(3) 



and the intensity (or the number of particles) 

A = / \ip(x,t)\ 2 dx 



(4) 



which, in this work, for convenience, we have set equal to one. Often we can "undo" this 
setting of N = 1 by an appropriate scaling of x, t and ip. 



III. COULOMB-TYPE INTERACTION 

First we consider a Coulomb-type interaction, i.e. the case of K(x,x') = —\x — x'\. 

Clearly the equation (J2J) is too hard to solve analytically. However, we expect it to have 
solitonic solutions. Note that both terms in (j3J), due to the — sign in K, are positive but 
their roles are very different; the first term tends to spread the field while the second one 
tends to focus it. Hence we expect the interplay of both terms to lead to a stable (stationary) 
solution which minimises the Hamiltonian H. 

What is this solution like? Clearly, to get some intuition for its shape we can try a 
gaussian of the form: 

^( x ) = ^l ex P (~Y x2 ) ' ( 5 ) 

If this approximation is relevant then the height of the soliton, defined here as the max of 
ipip*, in this Ansatz is proportional to a, which is a free parameter in (jSJ), while its width is 
proportional to 1/a. 
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Inserting (jSJ) into we find that 



^=k + ^- (6) 
The value of the parameter a which minimises the energy then is given by 

1/3 



a 



lit* 



(7) 



Note that increasing strength of the nonlocal interaction k, the height of the soliton (given 
by increases and the width (given by 1/a) decreases. For k = (the case of the "pure" 
i.e. linear Schrodinger equation), no localised structures appear, as expected. 

Note that, as expected, (JBJ) shows that the energy of the approximate solution consists of 
two terms which exhibit a different behaviour when one varies a. To predict the dependence 
of the minimal energy on the coupling k, we insert (J7J) into (jUJ) and find that 

H= 6 -^—-. (8) 

2(2?r) 1 /3 v ; 



A. Numerical results 

Next we have checked how good our predictions based on the gaussian approximation 
really are. To do this we have performed several numerical calculations using a fourth-order 
Runge-Kutta method for simulating the time evolution. Of course the x-variable has been 
discretised and we have typically used 500 lattice points. 
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FIG. 1: The value of if)ip* in dependence on x is shown for a typical soliton solution with k = 0.003. 

We have used several starting field configurations. Most of our simulations involved 
starting with a strongly localised structure, e.g. \ip(xo)\ 2 = 1, \if)(x ^ Xq)\ 2 = with 
Xq = 250. Then we took the energy out of the system by adding a small dissipative term 
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into (J2J having first taken most of the explicit time dependence of i/j(x,t) into account by 
writing ip(x, t) as 

ip{x,t) = exp(iEt) g(x,t), (9) 

where E can be calculated from (j2J), and really only evolving the equation for g(x,t). 

In consequence, the system evolved into a static equation for g(x,t) = g(x), i.e. a 
stationary solution field equation for ip(x,t). 

A typical soliton solution is shown in Fig|T](it corresponds to k = 0.003). We have found 
that the final stationary configurations were very well localised even for small values of k. 

Next we studied the dependence of the height and of the energy of the stationary solution 
as a function of the nonlocal coupling constant k. Our results are shown in FigJ21 




FIG. 2: The height of the soliton (^ip*) m ax as well as its energy for a nonlocal Coulomb interaction 
as a function of the nonlocal coupling k. 

Clearly, we have found that soliton solutions exist for all values of the nonlocal coupling 

K,. 

Fitting our numerical data for the height of the soliton, we see that {ipip*) m ax ~ 0.4732 • 
K o.37i3 for K e [q . L4 ] and ^*) max « 0.4372 • k - 3394 for k e [0 : 0.1]. For the energy, we 
have H = 0.7608 • k - 6390 for k e [0 : 1.4] and H = 0.7968 • k 06630 for k e [0 : 0.1]. Thus for 
small values of k our numerical results agree extremely well with the predictions based on 
the gaussian approximation - given by the expressions in ((7|) and (jHJ). 



IV. YUKAWA-TYPE INTERACTION 

Next we have looked at a Yukawa-type interaction, for which we have K (x, x') = 
^exp(— j\x — x'\). Note that 7 = corresponds to the case of the Coulomb-interaction 
studied above. 

This time it is even harder to get any analytical results. So, again we look at a gaussian 
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approximation. Thus we insert (jHJ) into (JSJ) and find that 
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where erf is the errorfunction. The value of this integral has to be evaluated numerically. 
We find that e.g. 1(1) ps 1.643, 1(2) ps 1.056, J(l/8) ps 2.851, i.e. we see that for fixed 7, 
the value of the integral increases with increasing a. The energy thus decreases with the 
increasing height of the approximate soliton. 

To be able to give arguments on how the height of the soliton depends on the coupling 
7, we expand the exponential for small 7 as 
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Neglecting 0(7 2 )-terms, we obtain for the energy 

1 



H=-a z - — + 
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(13) 



where the second term on the rhs is a constant piece that is not important for the arguments 
that follow. The value of a that minimises the energy fulfills the equation 



4 k 7^ 
a -=ol H = . 

V2T 2 



(14) 



Since the expressions for the roots of this equation are complicated, we evaluate the required 

( \ 1/3 

value of a by comparing it to the solution in the case when 7 = i.e. cto = I -7= ) . Thus 
we put a = ao + P7 and inserting this into (fUj) we find that 



a 



a - 



6 



7 for 7 small 



(15) 



This tells us that for fixed values of k, the height of the soliton decreases with 7. Moreover, 
in the linear approximation in 7, the correction to the height doesn't depend on the strength 
of the nonlocal coupling k. 



A. Numerical results 

Again, as in the Coulombic case, we have studied this problem numerically using the same 
procedures as before. In all studied cases we have found solitonic - like stationary solutions. 
Next we have looked at the properties of these solutions. We have fixed k and have tried to 
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FIG. 3: The height of the soliton {i^ip*) m ax for a non-local Yukawa interaction as a function of 7 
for k = 0.2 and k = 0.05, respectively. 

determine the height and the energy of the stationary soliton as a function of the parameter 
7. We have found that in all cases the soliton height decreases as one increases 7. In Fig|3] 
we present our results for k = 0.2 and k = 0.05, respectively. As predicted by our analytic 
argument (|15jh we observe that the height of the soliton decreases linearly with 7 and that 
the two curves are parallel to each other indicating the independence of the decrease on k. 

In Fig0J we show the energy of the solutions as a function of 7 for k = 0.2 and k = 0.05. 
Note that when we have added the constant piece k/ (27) to the energy (compare also (H2J)), 
in order to make the energy finite for 7 = 0, the energy also decreases. For 7 = 0, the 
Yukawa interaction reduces to the Coulomb interaction, while for 7 —>■ 00 we obtain the 
case of the "pure" Schrodinger equation. We would thus expect the height and the rescaled 
energy to tend to the values found in the Coulomb case for 7 — > and our numerical results 
confirm this. Moreover, for 7 — > 00, there is no localised structure and so the height and 
the energy, as expected, tend to zero. 

V. CONCLUSIONS 

We have studied a class of nonlinear Schrodinger models in which the nonlinear term 
involves nonlocal interactions based on Coulomb or Yukawa like potentials. We have shown 
that, like the more familiar local nonlinear Schrodinger model, these models also possess 
solitonic-like stationary solutions for all values of the coupling constant. All solutions in the 
Yukawa case are less localised than the corresponding Coulombic solutions. As expected, 
when one increases the strength of the potential the solitons become more localised. 

The solitonic fields look a little like gaussians; in fact an approximation of them by a 
gaussian is very reliable. We have calculated various properties of the solitons based on this 
approximation and have found them to be in a good agreement with what we have seen in 
our simulations of the full equations. 

In a similar gaussian approximation was successfully applied to the problems involving 
a local nonlinear term (but involving derivatives) in many dimensions. Would this be the 
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FIG. 4: The energy of the soliton for a non-local Yukawa interaction as a function of 7 for k = 0.2 
and k = 0.05, respectively. Note that we have subtracted the term K/27 from the energy (compare 
©). 

case when we look at nonlocal terms in two or even three spatial dimensions? This problem 
is currently being looked at. 
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